Dissipation in circuit quantum electrodynamics: 
lasing and cooling of a low-frequency oscillator 
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Superconducting qubits coupled to electric or nanomechanical resonators display effects previously 
studied in quantum electrodynamics (QED) and extensions thereof. Here we study a driven qubit 
coupled to a low- frequency tank circuit with particular emphasis on the role of dissipation. When 
the qubit is driven to perform Rabi oscillations, with Rabi frequency in resonance with the oscillator, 
the latter can be driven far from equilibrium. Blue detuned driving leads to a population inversion 
in the qubit and lasing behavior of the oscillator ("single-atom laser"). For red detuning the qubit 
cools the oscillator. Fhis behavior persists at the symmetry point where the qubit-oscillator coupling 
is quadratic and decoherence effects are minimized. Here the system realizes a "single-atom-two- 
photon laser". 



INTRODUCTION 



Recent experiments on quantum state engineering with 
superconducting circuits realized concepts originally in- 
troduced in the field of quantum optics, as well as 
extensions thereof, e.g., to the regime of strong cou- 
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Josephson qubits play the role of two-level atoms while 
electric or nanomechanical oscillators play the role of the 
quantized radiation field. In most QED or circuit QED 
experiments the atom or qubit transition frequency is 
near resonance with the oscillator. In contrast, in the ex- 
periments of Refs. with setup shown in Fig. QJi), the 
qubit is coupled to a slow LC oscillator with frequency 
(wt/27t ~ MHz) much lower than the qubit's level split- 
ting (AE/2irh ~ 10 GHz). The idea of this experiment is 
to drive the qubit to perform Rabi oscillations with Rabi 
frequency in resonance with the oscillator, ~ ll>t- In 
this situation the qubit should drive the oscillator and 
increase its oscillation amplitude. When the qubit driv- 
ing frequency is blue detuned, the driving creates a pop- 
ulation inversion of the qubit, and the system exhibits 
lasing behavior ( "single- atom laser"); for red detuning 
the qubit cools the oscillator A similar strategy for 
cooling of a nanomechanical resonator via a Cooper pair 
box qubit has been recently suggested in Ref. [20j. The 
analysis of the driven circuit QED system shows that 
these properties depend strongly on relaxation and deco- 
herence effects in the qubit. 

In the experiments of Ref. [l[ an enhancement of the 
oscillator due to the driving was observed. However, first 
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FIG. 1: a) In the setup of Ref. [IJ an externally driven three- 
junction flux qubit is coupled inductively to an LC oscillator, 
b) A charge qubit is coupled to a mechanical resonator. 



attempts to explain the effect theoretically did not re- 
solve several issues 2l|, |22[. In the experiments, in order 
to minimize decoherence effects, the Josephson flux qubit 
was biased near the flux degeneracy point. At this sym- 
metry point also the coupling to the oscillator is tuned 
to zero, and the enhancement should vanish. Uncon- 
trolled small deviations from the symmetry point might 
lead to the observed effect [jj, [5l| , but this explanation 
has not been supported by experiments. Here we explore 
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an alternative, namely that a quadratic coupling to the 
oscillator near the resonance condition w 2lut is re- 
sponsible for the observed enhancement. In the following 
we will consider both linear as well as quadratic coupling, 
which dominate away from the symmetry point and at 
this point, respectively. The second unresolved prob- 
lem is the magnitude of the effect. The experiments [l| 
showed an increase by a factor 4 — 5 in the amplitude, 
i.e., 16 — 25 in the number of oscillator quanta. The 
theory of Ref. [2l|, valid in the perturbative regime, pre- 
dicts a much weaker effect. We obtain a strong effect 
as follows [l8j]: due to a detuning of the qubit driving 
a population inversion is created at the Rabi frequency, 
and the system becomes a "single-atom laser" at the res- 
onance VLr ~ u)t, or a "single-atom-two-photon laser" for 
fl R w 2lot 2^, 24 1 . In both cases the lasing threshold is 
reached for realistic system parameters, and the number 
of quanta in the oscillator is increased considerably. 

A related situation, called "dressed-state la sing " , had 
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been studied before in quantum optics 
present scenario differs from that one in so far as the res- 
onator modes are coupled to the low-frequency Rabi os- 
cillations rather than to the high-frequency Mollow tran- 
sitions. The Rabi frequency can be readily tuned to reso- 
nance with the oscillator, which should facilitate reaching 
the lasing threshold and a proper lasing state. A similar 
idea has been explored in Ref. [27] in connection with 
coupling of atoms. 

In experiments with the same setup as shown in 
Fig. QJi) but in a different parameter regime the mech- 
anisms of Sisyphus cooling and amplification has re- 
cently been demonstrated [7j. Due to the resonant high- 
frequency driving of the qubit, depending on the detun- 
ing, the oscillator is either cooled or amplified with a ten- 
dency towards lasing. The Sisyphus mechanism is most 
efficient when the relaxation rate of the qubit is close to 
the oscillator's frequency. In contrast, in the present pa- 
per we concentrate on the "resolved sub-band" regime 
where the dissipative transition rates of the qubits are 
much lower than the oscillator's frequency. 

Also in situations where the qubit, e.g., a Josephson 
charge qubit, is coupled to a nano- mechanical oscillator 
(Fig.QJ>) it either cools or amplifies the oscillator. On one 
hand, this may constitute an important tool on the way 
to ground state cooling. On the other hand, this setup 
provides a realization of what is called a SASER [281 ]. 

Lasing and cooling of the oscillator has also been ob- 
served in a slightly different setup, when the ac-driven 
qubit is replaced by a driven superconducting single- 
electron transistor biased near the Josephson quasipar- 
ticle cycle [H 0, EH- When the SSET is coupled to 
a nanomcchanical or electric oscillator it can be used to 
either cool the oscillator 



32, 33, 34, \3 



& ISTJJ or to 

produce laser-like behavior. The latter has recently been 
observed in experiments 38 1. 



THE SYSTEM 



The Hamiltonian 



The systems to be considered are shown in Fig. [T] 
A qubit is coupled to an oscillator and driven to per- 
form Rabi oscillations. To be specific we first analyze 
the flux qubit coupled to an electric oscillator (Fig. QJl) 
with Hamiltonian 



H = --e($i c )a z - -Aa x -hn R0 cos(u> d t)(Tz 
+ huiTa a + ga z (a + a i ) . 



(1) 



The first two terms describe the qubit, with Pauli ma- 
trices a x . z operating in the flux basis of the qubit. The 
energy bias between the flux states e($!£ c ) is controlled 
by an external DC magnetic flux, and A is the tunneling 
amplitude between the basis states. The resulting level 
spacing AE = y/ e 2 + A 2 typically lies in the range of 
several GHz. The third term accounts for the driving of 
the qubit by an applied AC magnetic flux with ampli- 
tude CIro and frequency LOd- The last two terms describe 
the oscillator with frequency lot = 1/\/LtCt, which for 
the experiments of of Ref. [l| lies in the range of several 
10 MHz, as well as the qubit-oscillator interaction. Wc 
estimate the coupling constant g w MI p It,o to be of the 
order of 10 MHz. Here M is the mutual inductance, I p 
the magnitude of the persistent current in the qubit, and 
It.o = v IujJtI^Lt the amplitude of the vacuum fluctua- 
tion of the current in the LC oscillator. 

After transformation to the eigenbasis of the qubit, 
which is the natural basis for the description of the dis- 
sipation, the Hamiltonian reads 



H 
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AE a z — hfljio cos (u)dt) (sin £ a z — cos Q a x 



+ hu)Ta)a + g (sin ( a z — cos £ <? x ) ( a + o ) , (2) 



with tanC = e/A and AE = Ve 2 + A 2 . 

Because of the large difference of the energy scales 
between the qubit and the oscillator, AE Hojt, it 
is tempting, in the spirit of the usual rotating wave 
approximation (RWA), to drop the transverse coupling 
term — gcos(a x (a + a*) of Eq. ((2|). However, near 
the symmetry point (where sin £ = 0) the longitudinal 
coupling is weak. Therefore, we retain the transverse 
coupling, but transform it by employing a Schrieffer- 
Wolff transformation, Us — exp(iS), with generator 
S = (g/AE)cos( (a + at) a y , into a second-order lon- 
gitudinal coupling. On the other hand, since hiod ~ AE, 
we can drop within RWA the longitudinal driving term 
— hflfto cos (u/dt) sin C a z . The Hamiltonian then reads 

H = — — AE(T Z + hfljio cos (uJdt) cos £ a x + Hujt a* a 

g 2 
+ g sin C cr z (a + a f ) - — cos 2 Ca z (a + a)) . (3) 
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A further unitary transformation with Ur = 
exp (—iujdcr z t/2) brings the Hamiltonian to the ro- 
tating frame, H = U R HU R + ihil R U R . We obtain 

H = — hSuj <j z + — Ml no cos £ a x + Huit o) a 

+ gsin(a z (a + a 1 ") - ■^cos 2 Ccr z (a + a 1 ) . 



(4) 



Here 8ui = ujd — AE/h is the detuning. After diagonal- 
ization of the qubit terms we obtain 

1 



H = — Ml n cr z + hwT a) a 

+ g sin C [ sin (3 a z — cos /3 a x ] [a + a 
9 2 



cos 2 C [sin (3<j z - cos (3a x ] (a + a*) .(5) 



Here £Ir = y/ fl R0 cos 2 ( + 8lo 2 and tan/3 = 
Slo/(Qro cosC). 

Finally we employ a second RWA. While the first 
one dropped terms oscillating with frequencies of order 
AE/h, the second one assumes the Rabi frequency Qr 
and the oscillator frequency lot to be fast. In the inter- 
action representation with respect to the non-interacting 
Hamiltonian, H — {MlR/2)a z + hiOT a^a, we then obtain 

Hj = 5 i (aV.e-'^-^' 1 + /i.e.) 

+ g 2 (a) 2 a^e-^ R - 2 ^ t + /i.e.) 

+ gz (a) a + aa)) a z . (6) 

We kept both single-photon and two-photon interactions 
with gi = g sin £ cos (3 and g 2 — (g 2 / AE) cos 2 Ccos/3, 
although within RWA only one of them survives: the 
single-photon term for Qr ~ lot, or the two-photon 
term for Ur ~ 2u>t- The last term of © with 33 = 
— {g 2 / AE) cos 2 £sin/3 is the ac-Stark effect, causing a 
qubit state dependent frequency shift of the oscilla- 
tor [3§| . In what follows we will assume that the qubit is 
kept near the symmetry point, i.e., e <C A and cos£ ~ 1. 



Transition rates in the rotation frame 

The transformation to "dressed states" in the rotat- 
ing frame modifies the relaxation, excitation and deco- 
herence rates as compared to the standard results. To 
illustrate these effects and justify the treatment of the 
dissipation in latter sections we first consider a driven 
qubit (ignoring the coupling to the oscillator) coupled to 
a bath observable X B , 



H 
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AE a z + HD,rq cos (uidt) <j x 



(b x a x + byOy + b z a z ) X B + ffbath • (7) 



In the absence of driving, £Ir = 0, and for regular 
power spectra of the fluctuating bath observables we can 
proceed using Golden rule type arguments. The trans- 
verse noise is responsible for the relaxation and excitation 
with rates 

l -(X%(u = AE)) 



r 



r T = 



16-1 



Ah 2 

M 

m 2 



(X&u = -AE)) , 



(8) 



while longitudinal noise produces a pure dephasing with 
rate 



2h 2 



Sx(co = 0) 



(9) 



Here 6j_ = b x + ib y , and we introduced the ordered cor- 
relation function {X%{w)) = J dt e 4Wt (X B (t)X B {0)) , as 
well as the power spectrum, i.e., the symmetrized cor- 
relation function, S x (u) = {{X 2 b {lo)) + (X B (-u)))/2. 
The rates ([5]) and ^ also define the relaxation rate 
1/Ti = Ti = + Tj and the total dephasing rate 
I/T2 = T v = Ti/2 + r* which appear in the Bloch equa- 
tions. 

To account for the driving with frequency u>d it is con- 
venient to transform to the rotating frame via a unitary 
transformation Ur — exp (—iuid<J z t/2). Within RWA the 
transformed Hamiltonian reduces to 



H = ^h[Q R0 a x 



5uJ(7 z 



(10) 



2 [ b *<r* 



bj_e iWdt a- 



X f 



balli 



where b± = b x +ib y , and the detuning is Slo = LOd—AE/h. 
The RWA cannot be used in the second line of (fit)]) since 
the fluctuations X B contain potentially all frequencies, 
including those of order ±u>d which can compensate fast 
oscillations. 

Diagonalizing the first line of (fTU|) one obtains 



H = ^MIr<t z 



sin/3 , 

b z 

2 

(sin (3 - 



+ -Hbath 
COS/3 



(61 



1) 



,-iu>At + h± e iu d tj 

t , (sin/3-1) 



(T z Xf 



b± e 



iu) d t 



cos/3 



(j-i- Xr + h.c 



}■ 



(11) 



with Qr = \J Qrq + Slo 2 and tan/3 = Slo/^Irq. 

From here Golden-rule arguments yield the relaxation 
and excitation rates in the rotating frame [!(| 

- b 2 z cos 2 (3 



Ah 2 

\b±l 



{x B {Lu = n R )) 



(12) 



+ 



16ft 2 



(l-sin/3) 2 (X 2 B (LU = Lo d + nR)) 
(l + shi(3) 2 (X 2 B (LU = -LOd + nR)) 
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and 



b 2 z cos 2 /3 ! ^ 2 
4ft 2 

M 
M 

16ft 2 



(I 2 ( W = -Q a )) (13) 
(l-sin/?) 2 <.X-|(w = -w d -fi B )} 
(1 + sin/3) 2 <Xl(w = w d -n fl )) , 



as well as the "pure" dephasing rate 

- ^Isin 2 /? n ^ , cos 2 /? 2 

r,„ w — — — Sx(u = 0) + | &l| Sx(w = Ud) ■ 



2ft 2 



(14) 



We note the effect of the frequency mixing, and due to 
the diagonalization the effects of longitudinal and trans- 
verse noise on relaxation and decoherence get mixed. In 
addition we note that the rates also depend on the fluc- 
tuations' power spectrum at the Rabi frequency. 

For a sufficiently regular power spectrum of fluctua- 
tions with us as ±AE/h, we can ignore the effect of de- 
tuning and the small shifts by ±H, R as compared to the 
high frequency io d as AE/H. We further assume that 
£Ir <^ kgT /ft. In this case we find the simple relations 



1 t 
fx 



(1 + sin/3) 2 (1- sin/3) 2 1 2 
4 r l + 1 r T + g cos r « ' 

(1-sin^) 2 (l + sin/3) 2 1 2 
1 1 1 H 4 1 T + - cos /3 1 „ , 

sin 2 /?r; + ^(r i + r T ) , (15) 



where the rates in the lab frame are given by Eq. {5J and 
the new rate, T„ = | & 2 depends on the power 

spectrum at the Rabi frequency. 

To proceed we concentrate on the regime relevant for 
our system. Near the symmetry point of the qubit we 
have b z as 0. At low temperatures, k-sT <§; AE as hujd, 
we can neglect Tj as it is exponentially small. Thus we 
are left with 



1/T 



(IT sin/?)' 



cos 2 /3 



(16) 
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FIG. 2: Dressed states of a driven qubit near resonance. Here 
m is the number of photons of the driving field, which is 
assumed to be quantized. 



In a more careful analysis, not making use of the "white 
noise" approximation, we obtain for (3 = 



fx = (x 2 B (uj = uj d + n R )) 
f T (x 2 B (io = uj d -n R )) 



(18) 



For instance, for Ohmic noise and low bath temperature 
this reduces to Ti/T-[ ~ 1 + 2Q R /uJd- This corresponds 
to an effective temperature of order 2ftwd/fcB ~ 2AE/k-Q, 
which by assumption is high but finite. The infinite tem- 
perature threshold is crossed toward negative tempera- 
tures at weak blue detuning when the condition 



(l+sin/3) 2 
(1 - sin/3) 2 



2n R 



(19) 



where 



is satisfied. 

To illustrate how the population inversion is created 
for blue detuning we show in Fig. [5] the level structure, 
i.e., the formation of dressed states, of a near- resonantly 
driven qubit. For the purpose of this explanation the 
driving field is quantized. This level structure was de- 
scribed first by Mollow [26j. The picture also illustrates 
how for blue detuning a pure relaxation process, Tj = T , 
in the laboratory frame predominantly leads to an exci- 
tation process, T|, in the rotating frame. 



Tn = 



4ft 2 



(A >2 (^ = As/ft))«r i 



(17) 



The ratio of up- and down-transitions depends on the 
detuning and can be expressed by an effective tempera- 
ture. Right on resonance, where /3 = 0, we have = 
corresponding to infinite temperature or a classical drive. 
For "blue" detuning, /3 > 0, we find F-\ > Fj., i.e., neg- 
ative temperature. This leads to a population inversion 
of the qubit, which is the basis for the lasing behavior 
which will be described below. 



The Liouville equation in the rotating frame 



Within the approximation leading to the rates (|15| the 
Liouville equation governing the dynamics of the density 
matrix in the rotating frame can presented in a simple 
Lindblad form. Start in the lab frame the dissipation is 
accounted for by two damping terms, 



[H, p}+ Lq p + L R p . 



(20) 
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The qubit's dissipation is described by 

r. 



■>qP 



^- (2tj-p<7+ — pa+<7- — <J + (T-p) 



— (2(7 + p<7_ - pa-Vj, 

r* 

[o z pa z - p) , 



<T-<J + p) 



(21) 



with rates given by ([5]) and © ■ The resonator damping 
with strength parametrized by k, can be written as 41 1 



LrP = ^ (^ th + •"■) { 2a P at ~ atfl P _ pat a) 
+— N t h {2a) pa — aa) p — paa)} . 



(22) 



Here Nth = 1/ [exp(Swr/fcBT) — 1] is the thermal distri- 
bution function of photons in the resonator. 

The transformations to the rotating frame and RWA 
described above transform the Liouville equation to a 
new form, which in our approximation again has a Lind- 
blad form. In the interaction representation it is 



P = —h 



Hi,p 



+ L Q p + L R p , 



(23) 



where 



L QP = Y (2<7_po+ - po+cr_ - CT+o-_p) 



r 
1 



T 

- (2<7+ptT_ — pO-G+ — O-O+p) 

p) , 



(24) 



while the oscillator damping term is not affected by the 
transformation. Although at low temperatures in the lab 
frame the relaxation processes dominate (see Eq. I17[) , the 
transformation to the rotating frame introduces excita- 
tion and pure dephasing processes. 



THE SINGLE-QUBIT LASER 

In the following we will consider two resonance situ- 
ations, Qr ~ lot or ~ 2lot, when either the one- 
or the two-photon interactions dominate, and investigate 
the effects of blue or red detuning, Su> = cod — AE/h, of 
the qubit driving frequency. We also study the effects of 
detuning of the Rabi frequency flu relative to that of the 
oscillator. 



One-photon interaction 

When the Rabi frequency is in resonance with the os- 
cillator, flu w lot, the Hamiltonian ^ in RWA reduces 



to 

Hi = g x (a+ a- e -*(n*-<r)* + h.c\ 

+ g 3 (era + aa)) a z . (25) 

From here we can proceed in the frame of the standard 
semiclassical approach 4l|, |42j of laser physics with the 
following main steps: In the absences of fluctuations the 
system is described by Maxwell-Bloch equations for the 
classical variables a = (a), a* = (a'), s± = (<r±) and 
s z = (cr z ), which can be derived from the Hamiltonian 
(f2"B")) if all correlation functions are assumed to factorize. 
Next the qubit variables can be adiabatically eliminated 
as long as k, gi < T^r^, which leads to a closed equa- 
tion of motion for a. If we finally account for fluctuations, 
e.g., due to thermal noise in the resonator, a becomes a 



stochastic variable obeying a Langevin equation [4l| . 



C 



+ m • (26) 



Here C = 2g 2 , sf = -D / (l + |a 2 |/n ) is the sta- 
tionary value of the population difference between the 

qubit levels, and D Q = (f j — f -H /f i is the normal- 
ized difference between the rates with Pi = ff + Pj.. 
We further introduced the photon saturation number 
ri = t,pti/Ag\ and h = n (l + <5fi 2 /f 2 ), and the 
total dephasing rate = Ti/2 + T* . The detuning 
of the Rabi frequency enters in combination with a fre- 
quency renormalization, 5fl = flu — lot + ff3M 2 . The 
Langevin force due to thermal noise in the oscillator sat- 
isfies (£(*)£*(*')) = KN th 6{t - t') and (£(«)£(*')) = 0. 
Noise originating from the qubit can be neglected pro- 
vided the thermal noise is strong, reiVth ^> gf/T^. 



Two-photon interaction 

The two-photon effect dominates near the resonance 
condition flu w 2lot- In RWA the Hamiltonian reduces 
to 

Hj = g 2 (a t2 a_e-^- 2 ^* + /i.c.) 

+ 53 (a) a + aa)) u z . (27) 

The corresponding Langevin equation for the resonator 
variable reads 



a 



C 



T v + iSfl 



sf + 4ig 3 s? 



2 + ( 28 ) 



i.e., is of the same form as Eq. (f2T))) but with C = 4g 2 \a\ 2 
and s s z = —Do/ (l + (|a 2 |/no) 2 ). The photon saturation 
number is now given by no — (r^Ti/Ag 2 ) 1 / 2 : and no = 
no(l + 8fl 2 /V 2 )) 1 / 2 . Again represents thermal noise, 
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FIG. 3: Average number of photons in the resonator as func- 
tion of the driving detuning 8uo and amplitude f2.Ro . Peaks at 
Slu > correspond to lasing, while dips at 5u) < correspond 
to cooling. The inner curve corresponds to the one-photon 
resonance which exists only away from the symmetry point. 
Here we assumed e = 0.01A. The outer curve describes the 
two-photon resonance, which persists at e = 0. In domains of 
bistability the lowest value of n is plotted (leading to the sharp 
drops in both curves). We chose the following parameters for 
the qubit: A/2tt = 1 GHz, e = 0.01A, F /2tt = 125 kHz, the 
resonator: ujt/2tt — 6MHz, k/2tt = 0.34 kHz, and the cou- 
pling: g/2n — 3.3 MHz. The bath temperature is T = 10 mK. 



while noise arising from the qubit can be neglected if 
nN t h g\hjY^. The detuning of the Rabi frequency 
for two-photon interaction is given by Sfl = f2# — 2u>t + 
SsM 2 ■ 



Results obtained from the Langevin equation 

If one neglects the frequency shifts of the oscillator, 
i.e., for g% = 0, the Fokker-Planck equations correspond- 
ing to the Langevin equations ([26]) and (|28|) have exact 
analytic solutions 0. Also for g 3 ^ the Eqs. (j2l))) 
and (|28p written as a = —f(n)a/2 + £(t) can be trans- 
formed to equations for the average number of photons 
(|a| 2 ) = n in the form n — — (nRe[/(n)]) + nN t h- In 
the steady state, for n 3> 1 they can be approximated by 
fiRc[/(n)] = nN t h- The results of this analysis are shown 
in Fig. [3] To demonstrate both the one-photon and the 
two-photon effects we have assumed a small deviation 
from the symmetry point, e = 0.01A. The two-photon 
resonance (the outer one) persist even for e = 0, while 
the one-photon resonance (the inner one) vanishes there. 
We observe that the solution shows bistability bifurca- 
tions (see below). As a result we see in in Fig. [3] sharp 
drops of n for both resonances as only the lowest stable 
value is plotted . 

We can estimate the asymptotic solutions analytically. 
In the one-photon case, assuming n 3> ho, we obtain 




j-\ — i — i — i — i — =-. — i — i — i — i — i — i — i 

-15 -10 -5 5 10 15 

Detuning frequency (MHz) 



FIG. 4: Average number of photons n versus the detuning. 
The blue curves are obtained from the Langevin equations 
(|26[) and ()28|) . They show the bistability with the solid curve 
denoting stable solutions, while the dashed curve denotes the 
unstable solution. The red curve is obtained from a numerical 
solution of the master equation (|23[l . The driving amplitude 
is taken as £Iro/2it = 5MHz. The parameters of the qubit: 
A/2?r = 1 GHz, e = 0.01A, r /27r = 125 kHz, the resonator: 
uj t /2tt = 6MHz, k/2tt = 1.7 kHz, N th = 5, and the coupling: 
g/2-K = 3.3 MHz. 

from Eq.flU 

AT i ( _ - D o)Tl , . 

n ~ A^th H r ■ (29) 

2k 

This result holds independent of whether the second con- 
tribution due to the qubit is larger or smaller than the 
thermal number A'th as long as h 3> ho. In the two- 
photon case, assuming n 3> ho, we obtain from Eq. (|28p 

at i ( _ - D o)ri . . 

n — A th H ■ (30) 

K 

Solution of the master equation 

We also solved the full master equation (f2"3"]) numeri- 
cally, which provides access not only to the average num- 
ber of photons in the oscillator, n, but also to the whole 
distribution function P(n). To reach convergence with 
a limited number of photon basis states (n < 100) we 
assumed a low thermal number, A^h = 5 and a relatively 
high relaxation constant of the oscillator k/27t = 1.7 kHz. 
In Fig.|4]the solutions of the Langevin equations ([26]) and 
(|28jl and those of the master equation (|23|) is compared. 

In Fig. [5] the Fano factor F = ((n 2 ) - (n) 2 ) /(n) of 
the photon number distribution is presented. We observe 
two phenomena. First, in the regime of the lasing without 
bistability the Fano factor is reduced as compared to that 
of the thermal state. In the bistable regime it is increased 
due to the switching between the two solutions. 

In Fig. [S] the distribution function, P(n), for the num- 
ber of photons in the oscillator is plotted both for the 
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FIG. 5: The Fano factor (blue), and the average photon num- 
ber n (red). The parameters are as in Fig. [3] 
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FIG. 6: The distribution function, P(n), obtained by numer- 
ically solving the master equation (|23|l . Blue curve: cooling 
regime of the one-photon resonance with Qro — 2n x 5MHz 
and Slo = — 2n x 3.2MHz. Red curve: lasing regime of the two- 
photon resonance with Qro = 2n x 5MHz 8lo = 2n x 11.7MHz. 
We observe a peak in the P(n) distribution between n = 20 
and 30 as a result of the lasing behavior. Black curve: thermal 
distribution with jV t h = 5. The parameters are as in Fig. [3] 



cooling and enhancement regime. For comparison also 
the thermal (Bose-Einstein) distribution is plotted. 

The enhancement and cooling effects described here 
rely crucially on the transition rates, i.e., the dissipative 
effects as described by the Liouville equations. In or- 
der to demonstrate this dependence we show in Fig. 
the dependence of the average photon number on the 
qubit's relaxation rate at the one-photon resonance. We 
note a non-monotonic dependence on the qubit's relax- 
ation rate. Above the saturation threshold for n > no 
the pumping rate is limited by To, leading to a roughly 
linear growth of the photon number with increasing To, 
consistent with Eq. (|29|) . At the saturation threshold for 
n ~ no the effective coupling is determined by g\ and the 
photon number becomes insensitive to small variations of 





— numerical solution 
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FIG. 7: Average number of photons in the resonator as func- 
tion of the qubit's relaxation rate, To at the one-photon res- 
onance, Qr — lot for gz — and Nth = 5. The dark blue 
line shows the numerical solution of the master equation, the 
light blue solid line represents the solution of the Langevin 
equation, Eq. ([26}. The green and red dashed curves repre- 
sent respectively the saturation number no and the thermal 
photon number Ath. The parameters are as in Fig. [4] (except 
for To). 



Tq ■ Finally, for n < no , an increase of To predominantly 
increases the dephasing rate T v . As can be seen from Eq. 
(|26|) this destroys the coherent coupling between qubit 
and oscillator and the photon number decreases towards 
N t h- In Figure [7] we plot both the results of a numer- 
ical solution of the master equation and the Langevin 
approximation. We find good agreement between both 
(except for the bistability in the corresponding parame- 
ter regions). 



RESULTS AND DISCUSSION 

We summarize our main conclusions. Our results for 
the number of photons n are plotted in Fig. [3] as a func- 
tion of the detuning Slo of the driving frequency and driv- 
ing amplitude Qro- It exhibits sharp structures along two 
curves corresponding to the one- and two-photon reso- 
nance conditions, ^Ir — lot — 9an and Qr = 2lot — 93^. 
Blue detuning, Slo > 0, induces a strong population in- 
version of the qubit levels, which in resonance leads to 
one-qubit lasing. In experiments the effect can be mea- 
sured as a strong increase of the number of photons in the 
resonator above the thermal values. Red detuning pro- 
duces a one-qubit cooler with resulting photon numbers 
substantially below the thermal value. 

The bistability of the solution of the Langevin descrip- 
tion is illustrated in Fig. 2J In the range of bistability we 
expect a telegraph-like noise corresponding to the ran- 
dom switches between the two solutions. 

Potentially useful applications of the considered 
scheme are the lasing behavior and the creation of a 
highly non-thermal population of the oscillator as well as 
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the cooling. Within the accuracy of our approach we esti- 
mate that a population of order n = 1 can be reached for 
optimal detuning. A more detailed analysis is required 
to determine the precise cooling limit. 

So far we described an LC oscillator coupled to a 
flux qubit. But our analysis equally applies for a nano- 
mechanical resonator coupled capacitively to a Joseph- 
son charge qubit (see Fig. [lb). In this case a z stands 
for the charge of the qubit and both the coupling to the 
oscillator as well as the driving are capacitive, i.e., in- 
volve a z . To produce the capacitive coupling between 
the qubit and the oscillator, the latter could be metal- 
coated and charged by the voltage source V x . The dc 
component of the gate voltage V g puts the system near 
the charge degeneracy point where the dephasing due 
to the 1// charge noise is minimal. Rabi driving is in- 
duced by an ac component of V g . Realistic experimen- 
tal parameters are expected to be very similar to the 
ones used in the examples discussed above, except that a 
much higher quality factor of the resonator (~ 10 5 ) and 
a much higher number of quanta in the oscillator can be 
reached. This number will easily exceed the thermal one, 
thus a proper lasing state with Poisson statistics, appro- 
priately named SASER [28j |. is produced. One should 
then observe the usual line narrowing with line width 
given by «;iV t h/(4n) ~ K 2 A^ t h/Ti. Experimental obser- 
vation of this line-width narrowing would constitute a 
confirmation of the lasing/sasing. 
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